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ABSTRACT 

Results for the gravitational Wilson loop, in particular the area law for large loops in the strong 
coupling region, and the argument for an effective positive cosmological constant, discussed in a 
previous paper, are extended to other proposed theories of discrete quantum gravity in the strong 
coupling limit. We argue that the area law is a generic feature of almost all non-perturbative lattice 
formulations, for sufficiently strong gravitational coupling. The effects on gravitational Wilson loops 
of the inclusion of various types of light matter coupled to lattice quantum gravity are discussed as 
well. One finds that significant modifications to the area law can only arise from extremely light 
matter particles. The paper ends with some general comments on possible physically observable 
consequences. 



1 Introduction 



The identification of possibie observables is an important part of formulating a theory of quantum 
gravity. In general it is expected that these quantum observables will be represented by expectation 
values of operators which have physical interpretations in the context of a manifestly covariant 
formulation. In this paper, we focus on the gravitational analog of the Wilson loop [1,2], which 
provides physical information about the parallel transport of vectors, and therefore on the effective 
curvature, around large, near-planar loops. We will extend the analysis of earlier work [3, 4] to more 
general theories of discrete quantum gravity. A recent complementary discussion of the significance 
of physical observables in a quantum theory of gravity can be found for example in [5]. 

In classical gravity the parallel transport of a coordinate vector around a closed loop is described 
by a rotation, which is a given function of the affine connection along the space-time path. Then 
the total rotation matrix U(C) is given by the path-ordered (V) exponential of the integral of the 
affine connection via 



The gravitational Wilson loop then represents naturally a quantum average of a suitable trace (or 
contraction) of the above non-local operator, as described in detail in [3]. Its large distance (i.e. for 
loops whose size is very large compared to the lattice cutoff) behavior can be estimated, provided 
one makes some suitable assumptions about the short distance fluctuations of the underlying geom- 
etry, with the key assumption being the use of a Haar integration measure for the local rotations at 
strong coupling. A general result then emerges, which is that the Wilson loop generically exhibits 
an area law for sufficiently strong gravitational coupling (large G) and near-planar loops [3,4]. It 
should be noted here that in contrast to gauge theories, the Wilson loop in quantum gravity [6] 
does not provide useful information on the static potential, which is obtained instead from the 
correlation between particle world-lines [7,8]. Instrumental in deriving the results of [3] was the 
first order Regge lattice [9] formulation of gravity, discussed originally in [10]. 

Furthermore, from a semi-classical point of view, a vector's rotation around a large macroscopic 
loop is expected to be directly related, by Stoke's theorem, to some sort of average curvature 
enclosed by the loop. In this semi-classical picture one would write for the rotation matrix U 
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where A c v is the usual area bivector associated with the loop in question, 

Ag = \jdx»x v . (3) 

The use of semi-classical arguments in relating the above rotation matrix U(C) to the surface inte- 
gral of the Riemann tensor assumes (as is usual in the classical context) that the curvature is slowly 
varying on the scale of the very large loop. Then, in such a semi-classical description of the parallel 
transport process, one can re-express the connection in terms of a suitable coarse-grained, or semi- 
classical, Riemann tensor, and thus relate the quantum Wilson loop expectation value discussed 
previously to an observable large scale curvature. The latter is represented phenomenologically by 
the long distance, observed cosmological constant \ bs- 

It is important in this context to note, as an underlying theme, the close analogy between 
the Wilson loop in gravity and the one in gauge theories, both theories involving a connection 
as a fundamental entity. Furthermore, a lot is known about the behavior of the Wilson loop in 
non-abelian gauge theories at strong coupling, some of it from analytical estimates and some from 
large-scale numerical simulations. Let us recall that in non-Abelian gauge theories, the Wilson loop 
expectation value for a closed planar loop C is defined by [1] 

W{C) = < TvVexp{ig I A^x) dx^ > , (4) 

J c 

with Afj, = taA^ and the t a 's the group generators of SU (N) in the fundamental representation. In 
the pure gauge theory at strong coupling [1,2], it is easy to show that the leading contribution to 
the Wilson loop follows an area law for sufficiently large loops 

< W(C) > ~ eM-Ac/f) (5) 

A — » oo 

where Ac is the minimal area spanned by the planar loop C and £ the gauge field correlation 
length. Furthermore, it can be shown that the area law is fairly universal at strong coupling, in 
the sense that it is not too sensitive to specific short distance details of the 5C/(A r )-invariant lattice 
action. Indeed one expects the result of Eq. (5) to have universal validity in the lattice continuum 
limit, the latter being taken in the vicinity of the ultraviolet fixed point at gauge coupling g = 0. 

The fundamental renormalization group invariant quantity £ appearing in Eq. (5) represents 
the gauge field correlation length, defined for example from the exponential decay of connected 
Euclidean correlations of two infinitesimal loops separated by a distance |x|, 

G n (x) = < TrPexp{igf A^(x') dx'^( x ) TrPexp{ig<f A^x") dx""}(0) > c . (6) 
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Here the C e 's are two infinitesimal loops centered around x and respectively, suitably defined on 
the lattice as elementary square loops, and for which one has at sufficiently large separations 

Gu(x) ~ exp(-M/£) • (7) 

\x\ — > oo 

Thus the inverse of the correlation length £ is seen to correspond, via the Lehmann representation, 
to the lowest gauge invariant mass excitation in the gauge theory, the scalar glueball. 

Through the renormalization group £ is related to the f3- function of Yang- Mills theories, with 
£ the renormalization group invariant obtained from integrating the Callan-Symanzik /3-function, 

C\g) = const. A exp (- J" , (8) 

with A the ultraviolet cutoff, so that £ is then identified with the invariant gauge correlation length 
appearing in Eqs. (5) and (7). 

In an earlier paper [3], we adapted the gauge definition of the Wilson loop to the gravitational 
case, specifically to the case of lattice gravity, and in the context of the discretization scheme due 
to Regge [9]. On the lattice, with each neighboring pair of simplices s,s + 1 one can associate a 
Lorentz transformation U^ u (s,s + 1), which describes how a given vector V 1 transforms between 
the local coordinate systems in these two simplices. This transformation is directly related to the 
continuum path-ordered (P) exponential of the integral of the local affine connection, with the 
connection here having support only on the common interface between two simplices. The lattice 
action itself only contains contributions from infinitesimal loops, but more generally one might 
want to consider near-planar, but non-infinitesimal, closed loops C. Along this closed loop the 
overall rotation matrix will be given by 

U^(C) = [J] U s , s+1 ]\. (9) 
sec 

In analogy with the infinitesimal loop case, one would like to state that for the overall rotation 
matrix one has 

U" V (C) « \e 5 ^ B ^V , (10) 



V 



where B„v{C) is now an area bivector perpendicular to the loop and 5(C) the corresponding 
deficit angle, which will work only if the loop is close to planar so that can be taken to be 
approximately constant along the path C. By a near-planar loop around the point P, we mean one 
that is constructed by drawing outgoing geodesies on a plane through P. 

The matrix U^ U (C) in Eq. (9) then describes the parallel transport of a vector round the loop 
C. If that is true, then one can define an appropriate coordinate scalar by contracting the above 
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rotation matrix U(C) with an appropriate bivector, namely 

W(C) = u aP (C) U af3 (C) (11) 

where the bivector, oj a p{C), is intended as being representative of the overall geometric features of 
the loop (for example, it can be taken as an average of the hinge bivector uj a p{h) along the loop). 
Finally, in the quantum theory one is interested in the quantum average or vacuum expectation 
value of the above loop operator W(C), as in the gauge theory expression of Eq. (4). 

The next step is to relate the so defined, and computed, quantum average to physical observable 
properties of the manifold. Indeed for any continuum manifold one can define locally the parallel 
transport of a vector around a near-planar loop C. Then parallel transporting a vector around 
a closed loop represents a suitable operational way of detecting curvature locally. Thus a direct 
calculation of the vacuum expectation of the quantum Wilson loop provides a way of determining 
an effective curvature at large distance scales, even in the case where short distance fluctuations in 
the metric may be significant. 

For calculational convenience, the actual computation of the quantum gravitational Wilson loop 
in [3] was achieved by using a slight variant of Regge calculus, where the contribution to the action 
from the hinge h is given not by the original Regge expression 

S h = - k A h 5 h , (12) 

with k = l/8irG, but instead by the modified form 

S h = k -A h ti[(B h + e h) (V h - U^ 1 )] , (13) 

where Ah is the area of the triangular hinge where the curvature is located, Bh (called Uh in [3,4]) 
is a bivector orthogonal to the hinge, e is an arbitrary multiple of the unit matrix and the 
product of rotation matrices relating the coordinate frames in the 4-simplices around the hinge. 
The motivation for this second choice was that analytical calculations could then be performed 
more easily in the strong coupling regime, using methods analogous to the ones used successfully 
for gauge theories [1,2]. Indeed it can be shown [3] that this second action contribution is equal to 

S h = - k A h sm(8 h ) , (14) 

independently of the parameter e, where 5h is the deficit angle at the hinge. For small deficit 
angles one expects this to be a good approximation to the standard Regge action, and general 
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universality arguments would suggest that the lattice continuum be the same in the two theories. 
The expectation values of gravitational Wilson loops were then defined by either 



< W{C) > = < tr(C/i U 2 U n ) > 



(15) 



or 



< W(C) > = < tv[{B c + e h) Ux U 2 U n ] > , 



(16) 



where the C/jS are the rotation matrices along the path, and, in (ii), Be is a suitable average 
direction bivector for the loop C, which is assumed to be near-planar. The values of < W > in the 
strong coupling regime (i.e. for small k) can then be calculated for a number of loops, including 
some containing internal plaquettes. It was found that for large near planar loops around n hinges, 
to lowest non-trivial order (i.e. corresponding to a tiling of the interior of the loop by a minimal 
surface) , 



where a + (3 = n, and A is the average area of the plaquettes. Then using n = Ac /A, where Ac is 
the area of the loop, the area-dependent first factor can be written as 



where we have set £ = [A/\ log (A; ^4/16) |] 1//2 . Recall that for strong coupling, k — > 0, so £ is real, 
and that the quantity £ is in principle defined independently of the expectation value of the Wilson 
loop, through the correlation of suitable local invariant operators at a fixed geodesic distance. 
In the following we shall assume, in analogy to what is known to happen in non-abelian gauge 
theories, that even though the above form for the Wilson loop was derived in the extreme strong 
coupling limit, it will remain valid throughout the whole strong coupling phase and all the way up 
to the non-trivial ultraviolet fixed point, with the correlation length £ — > oo the only relevant and 
universal length scale in the vicinity of the fixed point. 

The next step was to interpret the result in semi-classical terms. By the use of Stokes's theorem, 
the parallel transport of a vector round a large loop depends on the exponential of a suitably- 
coarse-grained Riemann tensor over the loop. So by comparing linear terms in the expansion of 
this expression with the corresponding term in the expression of the area law, one can show [3] 
that the average curvature is of order l/£ 2 , at least in the strong coupling limit. Since the scaled 
cosmological constant is a measure of the intrinsic curvature of the vacuum, this also suggests that 
the cosmological constant is positive, and that the manifold is de Sitter at large distances. 




(17) 



exp[(Ac/A) log(*A/16)] = exp(-A c /£ 2 ) 



(18) 
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The question now arises as to whether these results are peculiar to the particular formulation 
of discrete gravity used. This led to a study of other proposed formulations, most of which were 
written down more than twenty years ago. In this work we will show that where it seems possible 
to define and calculate gravitational Wilson loops, the same area law emerges, and automatically 
implies a positive cosmological constant. 

Another key question we will address is whether these results are affected in any way by the 
presence of matter. After all the universe is not devoid of matter, and the pure gravity results 
should only be considered as a first order approximation to the full quantum theory (in a spirit 
similar to the quenched approximation in non-abelian gauge theories). This will be discussed here 
again in the context of the Regge formulation of discrete gravity used in [3] , using the methods of 
coupling matter to gravity reviewed for example in [4]. 

An outline of the paper is as follows. In Section 2, we describe formulations of Einstein gravity 
as a gauge theory on a flat background lattice, and in Section 3, the MacDowell-Mansouri descrip- 
tion of de Sitter gravity, as transcribed onto a flat background lattice by Smolin. More recent 
developments of discrete gravity, spin foam models, are discussed briefly in Section 4, and Section 
5 contains mention of other relevant theories. We then turn to the effect of matter couplings on 
the gravitational Wilson loop results, and Sections 6, 7, 8 and 9 contain systematic discussions of 
scalar matter, fermions, gauge fields and the lattice gravitino. Regarding these matter fields, the 
main conclusion is that the previous results are not affected, unless there are near massless spin 
1/2 and spin 3/2 particles (i.e. whose mass is comparable to the exceedingly small gravitational 
scale Section 10 consists of some conclusions. 

2 Gauge- theoretical treatment of Einstein gravity on a flat back- 
ground lattice 

We will first look at formulations of Einstein gravity as a gauge theory on a flat hypercubical 
background lattice, and in particular expand on the work of Mannion and Taylor [11] and of 
Kondo [12]. In these cases, the standard machinery for calculating Wilson loops in lattice gauge 
theories [2] can be taken over without too many modifications. Although such formulations were 
not the first chronologically of those we consider in this paper, we treat them first because they are, 
in many respects, the simplest. The idea is to write Einstein gravity in four dimensions, without 
cosmological constant, on a flat background lattice, treating it as a gauge theory with gauge group 
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SL(2, C), and relating it to the Einstein-Cartan formalism. In fact, for simplicity, we shall consider 
a Euclidean version, replacing SL(2,C) by SO(4). 

In the following nearest neighbor sites are labelled by n and n + fx, and their frames are related 

by 

U^n) = exp(iA„(n)) = U-^n + fi)' 1 , (19) 

where 

A^n) = l ia Af(n)S ab , (20) 
with a the lattice spacing and S a b the 0(4) generators, represented by the 4x4 matrices 

Sab = \ ha, Jb] , (21) 

with the Euclidean gamma matrices, 7 a satisfying 

{7a, 7b} = 2<U, 7a = 7a, a=l,...,4. (22) 

The curvature round an elementary plaquette spanned by the fi and v directions is given as usual 
by 

U^(n) = U^(n)U u (n + n)U^{n + n + !/)[/_„(« + 1/) = U^(n)U u (n + /x)^(n + i/)- 1 */^)" 1 , (23) 
and it can be shown that in the limit of small lattice spacing, 

U^uin) pa exp(ia 2 i?^), (24) 

where 

V = ^-^ + i[^,^] . (25) 

One notices that the usual lattice gauge theory type action, consisting of sums of U^ v terms, would 
give an R ilv W v term in the limit of small a, so terms involving the vierbein e^(n) and the matrix 
75 = 71727374 have to be introduced. One defines 

S = Tlh E ^ ^E x (n)U^{n)E a {n)] (26) 

n,fi,u,X,p 

where E^{n) = oe^7„ and k is the Planck length in suitable units. It can then be shown that 

S =^2 E e^e abcd Rf u (n)e c x (n)ej(n)+0(a 6 ), (27) 

n,n,u,X,p 

which is the Einstein action in first order form [13]. Furthermore by construction the action is 
invariant under local 0(4) rotations. For reasons which will become apparent, we shall consider a 
symmetrized form of the action: for each plaquette, rather than having the E a ^E\ term inserted 
only at the base point, we shall consider the average of its insertion at all vertices of the plaquette. 
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Figure 1. Illustration of the gravitational analog of the Wilson loop. A vector is parallel-transported along 
the larger outer loop. The enclosed minimal surface is tiled with parallel transport polygons, here chosen 
to be triangles for illustrative purposes. For each link of the dual lattice, the elementary parallel transport 
matrices U(s, s') are represented by arrows. 

In the following the partition function is defined by the usual path integral expression 

Z = J[dA][dE] exp(-S), (28) 
where [dA] = Jln^dnU^n), [dE] = Y\ n ^\dE\{n)^ and dnU is the Haar measure on 50(4). 
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Figure 2. A parallel transport loop, spanned by the 7 and S directions, with four oriented links on the 
boundary. The parallel transport matrices U along the links, represented here by arrows, appear in pairs 
and are sequentially integrated over using the uniform measure. 



Our interest here is in the definition and evaluation of Wilson loops in the strong coupling expansion. 
The authors of Ref. [11] define the loop around one plaquette, spanned by the 7 and 5 directions 
(see Fig. 2), by 

W = J] e 5 ^ Tr[E K {n)U s ^n)E^n)} , (29) 
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and so 

<W>= ^ j [dA][dE] [] e 57< Ti[E K (n)U 5l (n)E ( (n)} exp(-,S) . (30) 

They go on to show that in the strong coupling expansion, the dominant term is proportional to 

<W> = J [dE] e 5 ^ e s ^ e ahst e s K e\ e\ e b a , (31) 

where there is no sum over 7 and 5. Now suppose that 7 = 1, 5 = 2. Then the sum over A and a 
leads to 

e abst e 3 e 4 e 3 e 4 > 

(32) 

which is zero on symmetry grounds. Therefore their definition needs some modification, or one 
has to go to higher orders in the strong coupling expansion. In the latter case, it is possible to 
get a nonzero contribution by going to order l//c 6 , but here we concentrate on the first possibility. 
Omitting the Es from W also gives zero for < W >, so the modification we make is to insert a 75 
into W. The lowest order contribution is then 

- y-Ut J[dA\[dE\Y,e 5 ^T,[ lb E K {n)U 5l {n)E^n)\ £ Tr[ 75 E x (n') U^(n') E c {n')\ 

k,£ n' ,fj,,f,X,p 

[dA] [dE] ]T e*"* Tr[ 75 7s U s (n) U^n + 5) U s (n + 7)- 1 U^n)' 1 7t ] 
x £ e^ ACT Tr[ 75 7a C/ M (n') + /i) U^n' + ^ ^(n')" 1 76 ] <(n) e* (n) e£(n') e^(n') . 

(33) 

The integration over the As is equivalent to the integration over U's in SO(4) with the Haar 
measure: 

J d H U Un U H l = \ 5u S jk , (34) 

and we obtain 

[dE] ]T e s ^ e s ^ Tr[ 7d757c7f) 7 57a ] e%(n) e b (n) e c K (n) ef (n) . (35) 
Now we compute 

Tr[ 7d757c76757a ] = 4 (5 ab S cd - 5 ac 5 bd - 5 ad S bc ), (36) 

and, using 

g aX = \ a 2 Tr( 7a7fc ) e a x e b a = a 2 5 ab e a x e b a , (37) 
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we obtain 

16k 2 „ 

K§ Act 

Suppose that 7 = 1,5 = 2, then the sum over the indices in the e's and g's gives 



^ J [dE] J2 ^' Xa (5Act<7< - 9x^. ~ 9\K9t*) ■ (38) 



4(534 - 533 544) • (39) 

We expand the metrics in terms of the vierbeins and define the measure of integration to include 
a damping factor (Aa 2 /7r) 8 exp[— Aa 2 Sh. M (e^) 2 ] at each point, with ReX > [14], obtaining 

< 40 > 

(Note that we are ignoring a possible factor of the determinant of the vierbein in the measure.) 



U 2 U 2 



A >' 



u, 



D 



A V 



Figure 3. A vertex where various parallel transport matrices enter and leave, and where there are insertions 
on their paths. 

Before considering larger loops, let us obtain an algorithm which simplifies the calculations con- 
siderably. Consider a vertex with the matrices A, B, C, D attached to it, and [/-matrices attached 
to the lines entering and leaving the vertex, as shown in Fig. 3. Integration over the Us of the 
expression 

(UiU A bc (U 2 ) cd (U^) ef B fg (U 3 ) gh (Us% C jk (Ui) kl (U^) mn D no {U^) op (41) 

gives 

1 Tr(ABCD) 5 de S hi S lm S pa . (42) 
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We see that the effect of the integration is to give a factor of \ for each U, and to give the trace 
of the product of factors at each vertex. For a vertex with no insertion, we obtain the trace of the 
identity matrix, 4, and for one insertion of E^E, the value is zero since it is traceless. For two 
insertions of E^sE, we obtain 12/A 2 , where the integration over the Es has been done. Recall that 
there is also a factor of — 1/16k 2 for each plaquette, corresponding to the relevant terms in the 
expansion of the exponential of minus the action. This means that within our loop, if it is to have 
a nonzero value, every vertex must have either no Ej^E factors or two of them. This is why we 
took the average of insertions at all vertices of the plaquettes in the action; it would be impossible 
to get non-zero contribution from the internal plaquettes otherwise. 

Before proceeding with the calculations, let us mention an alternative to the procedure of 
averaging the contribution of the action from a plaquette over its vertices, a possibility, similar to 
the procedure in [3]. If we replace 75 by 75 + 6/4 for some arbitrary parameter e, the continuum limit 
of the action acquires a term proportional to g^ Ao ' R^ v \ a , which is zero because of the symmetries 
of the Riemann tensor, so the action is unaltered. However, the value of the Wilson loop is still 
zero, not because of the traces but because of factors of the Kronecker delta, which give zero on 
symmetry grounds. 
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Figure 4. A parallel transport loop around two plaquettes, with insertions of -E75-E on the loop shown by 
large dots. 

For a loop around two plaquettes, we find that we obtain a nonzero value only if the Ej§E is 
inserted at the place where the loop meets the second plaquette (see Fig. 4). The value of < W > 
is then 

1 (J)?) 2 ■ < 43 > 
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Figure 5. A larger parallel transport loop with twelve oriented links on the boundary. As before, the parallel 
transport matrices along the links appear in pairs and are sequentially integrated over using the uniform 
measure. The new ingredient in this configuration is an elementary loop at the center not touching the 
boundary. As in Fig. 4, the insertions of Ej^E on the loop are shown by large dots. 

For a loop around many plaquettes, we choose to insert the factors of Ej^E in the loop wherever 
the loop meets a new plaquette. For a loop with internal plaquettes, there has to be an even number 
of internal plaquettes as the insertions need to be paired between them. (For example, see the loop 
around nine plaquettes in Fig. 5; there is no way the insertions on the one internal plaquette can 
give a nonzero value.) This means that we obtain nonzero values only for Wilson loops surrounding 
an even number of plaquettes; the simplest case, with twelve plaquettes, is shown in Fig. 6. There 
are two ways of getting nonzero values from the internal plaquettes, corresponding to pairings of 
the insertions at the two vertices they have in common, which gives a factor of 2 in the answer, 
which, when integrated over the vierbeins, is 

— f— V 2 (44) 

Larger loops can then be treated in a similar way. From the results obtained so far, we deduce 
that, as the authors of [11] claimed, there is indeed an area law for large Wilson loops. The 
physical interpretation is of course very different, as discussed in the introduction, and later in the 
conclusions. 
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Figure 6. A Wilson loop around twelve plaquettes, of which two are internal, with the insertions of E^E 
on the loop shown by large dots as before. 

We now consider briefly the work of Kondo [12]. His basic formalism is very similar to that 
of [11], except that rather than introducing the vierbeins into the action directly, he introduces 
exponentials of them, with the action 

S = -4^2 E ^ Tr[ l5 U^(n)H x (n)H a (n)] , (45) 

n,ii,u,X,p 

where 

Hfj,(n) = exp[ioe^(n)7 a ] . (46) 

This has the consequence that the action is bounded. (The minus sign, which appears different 
from the sign in the formalism of [11], is because of the different relative position of the 75 factor.) 
In practice, in calculations, it is impossible to work out traces without expanding the exponentials 
and retaining the lowest order terms in the lattice spacing, so the formalism reduces to that of 
[11] in this respect, and the same values are obtained for the Wilson loops. (We have checked that 
the lowest order contribution comes from the product of the linear terms in the expansions of the 
exponentials.) However, Kondo also aims to set up a formalism which has reflection positivity, so 
his action contains sums over reflections, and if this full action is used, it is very complicated to 
evaluate Wilson loops. 

Note that this method of averaging the action contribution of each plaquette over the vertices 
of the plaquette needs to be used here, and could also be used in [3], eliminating the necessity for 
introducing the parameter e. 
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3 Lattice formulation of MacDowell-Mansouri gravity 



An earlier version of lattice gravity was given by Smolin [15], who transcribed the MacDowell and 
Mansouri [16] formulation of general relativity onto a flat background lattice. MacDowell and 
Mansouri built a gauge theory by defining ten (antisymmetric) gauge potentials by 

Af = rf, 4 a = y<, (47) 

where tt/f and are the usual gravitational connection and vierbein, and I is a lattice spacing. 
The curvature and torsion are defined in terms of the gauge potentials, and the action is of the 
form 

S = Jd 4 x e abcd Rf v , (48) 

where R°^ v is the Riemann tensor for 0(3, 2) or 0(4, 1). This can be shown to be equivalent, after 
multiplication by =Fl/32/ 2 /re 2 , with re the bare Planck length, to 

/ 2 1 9 



(iv (XT r>Oab r>Ocd , _J_ R _ f 



(49) 



where i2°" 6 is the usual Riemann curvature tensor. The first term is a topological invariant, the 
Gauss-Bonnet term, while the second and third are obviously the Einstein term and the cosmological 
constant term respectively, with a scaled cosmological constant A = ±2/re 2 Z 2 . Note that in this 
formulation the relative coefficients of various action contributions are fixed in terms of the bare 
parameter re and I. 

Then the starting point in [15] is the continuum action 

S = T l f d 4 x R% D e ABCD5 , (50) 

where is the curvature associated with an 0(4, 1) (minus sign) or 0(3,2) (plus sign) gauge 
connection, 6abcd5 is the totally antisymmetric 5-tensor and g = k/1 a dimensionless coupling 
constant. The parallel transport operators along the links of the lattice are defined by 

f-n+(i 



C/ M (n) = V exp 



rn+(i 

\g\ dx?Af{x)T AB 

J n 



(51) 



where the T AB are matrix representations of the relevant Lie algebra. Then the curvature around 
a plaquette on a hypercubic lattice, U^ v {n), is identical to the definition of Mannion and Taylor 
[11] [ Eq. (23) ], and this is related to the curvature by 

\ [U^(n)h = a 2 gR A ^ (T AB ) l3 + 0(a 3 ). (52) 
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The continuum action is then transcribed onto the lattice as 

S = ^E ^ z ijkl5 P^(n)]ij [U P v{n)] kl e ABC D5 ■ (53) 

9 n 

It involves a sum over contributions from perpendicular plaquettes at each lattice vertex, in analogy 
to the construction of the FF term in non-abelian gauge theories. In order to maintain the discrete 
symmetries of the lattice (reflections and rotations through multiples of it), this is extended to a 
sum over all orientations of the dual plaquettes 

S = feEE z ijkl5 KM)h [E#(n)]« e A BCD 5 . (54) 

g n 0,0> 

The partition function is then given by 

Z = J [dU] exp(iS), (55) 

where we take [dU] to be the normalized Haar measure. We restrict the integration to 0(5), rather 
than considering also 0(3, 2) and 0(4,1) as in [15], since for the non-compact groups one has 
to define the measure by dividing through by the (infinite) volume of the gauge group. For the 
five-dimensional representations used, the relevant integrals are: 

j[d H U) = l , (56) 

J [d H U] [U^nfiij = , (57) 

J [d H U] [U^n)}^ [U v {n )] kl = - 5 a 5 jk 5 nn > 5^ . (58) 

The structure of the action, based on pairs of dual plaquettes, means that the calculations are 
somewhat different from the case in [11]. In particular, since we want eventually to evaluate 
Wilson loops for planar surfaces, we can take as our basic building block a combination of two 
pairs of dual plaquettes, put together so that one plaquette from each pair lies adjacent to the 
other in the plane or they meet at one point, and the other two are joined back-to-back (see Figs. 7 
and 8). We then calculate the contribution from this configuration in both cases, when integration 
over the Us on the back-to-back faces is performed. 
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A 




Figure 7. Two pairs of dual plaquettes joined together, with the ones in the (/i, ^)-plane lying side-by-side, 
and the ones in the (/?, (r)-plane back-to back. 



Figure 8. Two pairs of dual plaquettes joined together, with the ones in the (/i, i^)-plane sharing only the 
vertex n, and the ones in the (p, <r)-plane back-to back. 



The quantity to evaluate in the first case is 

S = /M m 2 ^^''' 5 [P (1 W^(" + ^ 1 (n + ")0; 1 a 

x [U p (n)U a (n + p)Uj x (n + a)U- l {n)] k i p v {n)U^\n + u- n)U~ x (n - ^U^n - n)]^, 
x[i7 p (n)i7 CT (n + p)i7- 1 (n + ( T)C/- 1 (n)] fcT , (59) 

(with no summation over fi, v). Integration over the U p s and U a s gives 

(^) 2 (e^n 2 £ m5 e jj ' m [^(n)^(n + / x)C/- 1 (n + I ,)^- 1 (n + z,-^)C/r 1 (n- / x)C/ M (n-^)]^. (60) 



Now 



e ijkl5 e jj'kl5 = 2 ,M' vij vjj' 



(61) 
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where 



(62) 



so the final contribution, including a factor of 4 from the summation over p and a, is 




x [C/- 1 ^ - |u) U- l {n - n-v) U^n - fi - v) U u (n - u)] iff . (64) 



We now define a Wilson loop as the product of the U factors around the given path, with no extra 
factors in this case, and we calculate its expectation value as usual: 



As explained in [15], calculations are done in this formalism on the assumption that one can ignore 
the zero-torsion constraint; the basis for this is that the torsion is suppressed by a factor of j, where 
I is large. As a result, one only needs to integrate over the C/'s, and there is no need to integrate 
over the vierbeins in this formalism. Note that because of the structure of the basic building blocks, 
we can define Wilson loops only around paths which contain an even number of plaquettes. The 
simplest of these is shown in Fig. 4, and the area can be tiled by only one of the two possible 
building blocks, giving the value (l/(16# 2 ) 2 ) (192/125). 

The next most simple cases are shown in Fig. 9. The first of these can be tiled in four possible 
ways with the first of the building blocks, giving (l/(16g 2 ) 4 ) (2 12 3 2 /5 7 ), while in the second, which 
can be tiled in eight ways with the first building block and in one way with the second, the final 
contribution is (l/(16g 2 ) 4 ) (2 8 321/5 7 ). For the simplest configuration with internal plaquettes, a 
loop surrounding 12 plaquettes (see Fig. 6), there are many (1072) different ways of tiling it, 
so we need to add the contributions from all the different ways. The tiling shown in Fig. 6 
gives (l/(16g 2 ) 12 ) (2 26 3 4 /5 21 ) , and then combining this with the other contributions, we obtain 
(1/(16# 2 ) 12 ) (2 30 3 4 9481/5 23 ) . Notice the dependence on 1/g 2 in the various cases evaluated. Again, 
larger loops can then be treated in a similar way although the calculations become increasingly 
tedious. This indicates the usual area law for the gravitational Wilson loop. We note here that 
the authors of Ref. [17] have performed numerical simulations using the action from [15], but they 
restricted themselves to the gauge group 0(5), and to the weak coupling regime. 




(65) 
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Figure 9. Two arrangements for Wilson loops around four plaquettes. 



A formalism related to that of Smolin is described by Das, Kaku and Townsend [18]. They 
transcribe West's de Sitter invariant formulation of Einstein gravity onto the lattice, obtaining an 
action with plaquette contribution proportional to the square root of the trace of a square involving 
Smolin's action. They showed that their theory agrees with the one in [15] in the lattice continuum 
limit. The square root in the action makes it almost impossible to do any general analytical 
calculations. 

To put the results of Wilson loop calculations in this and the previous section, together with the 
results of [3], into context, it is interesting to make comparisons by relating the coupling constants 
to that of the continuum action. The k of Mannion and Taylor [11] and Kondo [12], and the g of 
Smolin [15] are related to the k = 1/8itG of [4] by 

1 k 1 _ l 2 k 

^ ~ 2^' ^ - 32 ' ( 66 ) 
Making a further normalization of the constants involved by equating the results for the smallest 
loop results, the answers of [11] and [12] agree with those of [3] until the loops contain internal 
plaquettes, and then, for example, the twelve-plaquette results differ by a factor of A 2 /6. The 
results of [15] are of the same order of magnitude as those of [3]. 
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4 Spin foam models 



Spin foam models grew out of a combination of ideas from the Ponzano-Regge model of three- 
dimensional discrete Lorentzian quantum gravity, and from loop quantum gravity. In loop quan- 
tization, the fundamental excitations are loops created by Wilson loop operators analogous to the 
ones used in gauge theories [19], and one assumes that states can be written as power series in 
spatial Wilson loops of the connection [20]. What does this intimate connection between Wilson 
loops and spin foam models mean in the context of this paper? 

In the three-dimensional formulation of Turaev and Viro [21], which is a regularized version 
of the Ponzano-Regge model, it has been shown [22] that the graph invariant defined by Turaev 
[23] coincides, in the semi-classical limit, with the expectation value of a Wilson loop. This is 
a consequence of the asymptotic behavior of 6j-symbols, with certain arguments fixed, involving 
rotation matrices which combine to give parallel transport operators along the graph. The extension 
of this result to graph invariants in discrete four-manifolds has not been made (as far as we know) 
and it is not clear anyway whether an area law could be obtained for large loops since the concept 
of a planar loop is not well-defined. 

One way of obtaining a spin foam model is from BF theory [24]. In four dimensions, repre- 
sentation labels are assigned to triangles and group elements to sections of the dual loop around 
each triangular hinge. The integral of the group elements around the dual loop gives the holonomy, 
which is a measure of the curvature, F. Thus an evaluation of Wilson loops is a basic ingredient 
in calculating the action, which is then conventionally expressed in terms of sums over amplitudes 
for the vertices, edges and faces of the spin foams. Alternatively, in group field theories, the action 
involves the integral over products of functions of the group variables, corresponding to a kinetic 
term and an interaction term. Here the evaluation of Wilson loops is somewhat similar to the way 
matter is inserted; certain edges are picked out (to form the loop) and are then treated differently 
in the summation process [25]. 

The authors of Ref. [26] have shown that there is an exact duality transformation mapping 
the strong coupling regime of non-Abelian gauge theory to the weak coupling regime of a system 
of spin foams defined on the lattice. They obtain an expression for the expectation value of 
a non-Abelian Wilson loop (or spin network) in terms of integrals of expressions involving finite- 
dimensional unitary representations, intertwiners and characters of the gauge group, together with a 
gauge constraint factor for each lattice point. The integrals are done explicitly, leaving complicated 
products and sums over intertwiners, projectors and the character decomposition of the exponential 
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of the action. Their calculation is very general, and to evaluate a Wilson loop in the usual sense, 
considerable simplification can be made. The links which form the Wilson loop can all be labelled 
with the same representation, and, as the loop has no multivalent vertices, the intertwiners all 
become trivial. Even so, the calculation is very complicated for a general gauge group, 

To illustrate the ideas behind the work of these authors, we will describe the corresponding 
calculations in lower dimensions and with gauge group SU(2) [27, 28]. We shall summarize the 
description in [28] . The partition function for gauge theory on a cubic lattice is written as usual as 
an integral over link variables Ui, with the action being a sum over plaquettes contributions 

z = p J U du i -p { 2T7I J ' (67) 

with (3 being the dimensionless inverse coupling. The matrix U p \ is the standard product of four 
link matrices U\ around the plaquette. The idea of the duality transformation is to make a Fourier 
transform in the plaquette variables, by first inserting unity for each plaquette into the partition 
function, in the form 

1 = I] / dU vi S ( U Pli UiU 2 U 3 U 4 ), (68) 
Pi 

where Lq...4 are the link variables around the plaquette. The (5-function can be realized by products 
of Wigner D-functions 

8{U, V)= J2 ( 2J + !) D Lm, ^) < m . (Y) . (69) 

7=0,^,1,... 

The unity is then inserted into the partition function in the form 

i = n / du vi E( 2j + i )^ im2 ( [/t ) 
P i j j 

^Di^U^Di^^Di^U^Di^U,) . (70) 

The integration over the plaquette matrices is performed using 

/ dU pl exp ^ Epl ( ^ + C ' C ° ) D J mim2 (U^) = | 5 mi , m2 h(P)TjW), (71) 

where Tj((3) = hj+iifi) / h{fl) [2,29] is the "Fourier transform" of the Wilson action and the /„ 
are modified Bessel functions. The partition function is then 



| AW) 



no. of plaquettes _ 

j P pi 



IJ JdUiDk^WD^WD^WD^iUt). (72) 



x 

links 1 
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In two dimensions, each link is shared by two plaquettes and the integration over Ds gives Kronecker 
deltas, whereas in three dimensions, each link is shared by four plaquettes and the integration over 
L>s gives 6j-symbols as in the Ponzano-Regge model. To compute the expectation value of a Wilson 
loop in representation j s , a factor of D^ S (U) must be inserted for each link on the loop. In two 
dimensions, use of the asymptotics of Tj{(3) leads to the area law at large (3 (strong coupling) 
[28]. In three dimensions, the extra Ds along the link give rise to 9j-symbols, and the asymptotic 
behavior of the Wilson loop has not been calculated explicitly. 

The formulation of spin foam models which seems the most tractable for the calculation of 
Wilson loops is the one of Ref. [30]. (Their expressions are essentially identical to those written 
down earlier by Caselle, D'Adda and Magnea [10,3] (see also [31]). In the absence of a boundary, 
their action can be written as (see Eq. (13) ) 

S = Y, Tr[BfW f (t)} , (73) 
/ 

where the sum is over triangular hinges, /, Uf(t) is the product of rotation matrices linking the 
coordinate frames of the tetrahedra and four-simplices around the hinge and Bj(t) is a bivector 
for the hinge, defined as the dual of £/(t). This in turn is the integral over the triangle / of the 
two-form S(t) = e(i) A e(t), formed from the vierbein in tetrahedron t. The action is independent 
of which tetrahedron is regarded as the initial one in the path around the hinge. There is a slight 
subtlety in the definition of Uf(t), as the basic rotation variables are taken to be Vt v , which relates 
the frame in tetrahedron t to that in 4-simplex v, of which t is a face, which is crossed in the path 
around hinge/. Then 

U f (t) = V tvl V Vltl ...V Vnt . (74) 

The action is sufficiently similar to that used by us in an earlier paper [3], that we may take over 
the formalism for calculating Wilson loops from there. The integration over the Vs, which are 
elements of SO (4) in the Euclidean case, proceeds exactly as in [3], and the same problem arises 
with the unmodified action of [30], as the bivector B is traceless. Therefore the definition of the 
action and of the gravitational Wilson loop has to be modified by an addition of e/4, as in [3], 
which again does not affect the value of the action. The results obtained are equivalent to those 
in our earlier paper, which indicates that the area law also holds for this formulation of spin foam 
models. 
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5 Other discrete models of quantum gravity 



We now consider very briefly various other approaches to discrete quantum gravity and the pos- 
sibility of evaluating the expectation values of gravitational Wilson loops in them. Kaku [32] has 
proposed a lattice version of conformal gravity, with action 



where P^(n) gives the curvature round a plaquette and is related to the Us in our previous 
equations, with U^(n) given in terms of the 0(4, 2) generators. The strong coupling expansion of 
the partition function is given by 



where the last term is included to impose the zero-torsion constraint. The analytic calculation 
of Wilson loops is complicated considerably by the presence of this constraint. If it is ignored, 
and Wilson loops defined as a product of Us round the loop as usual, then comparison with other 
calculations suggests that an area law will be obtained. (The calculations are very similar to those 
of Ref. [15] if one assumes a form for the 0(4,2) integrals as in his paper. The 75s disappear in 
the process of evaluating the basic building blocks.) 

Rather than considering conformal gravity, Tomboulis [33] has formulated a lattice version of 
the general higher derivative gravitational action in order to prove unitarity. He uses the gauge 
group 0(4) and considers vierbeins coupled as "additional matter fields", as in Mannion and Taylor 
[11] and Kondo [12], together with further auxiliary fields. After including reflections in order to 
preserve discrete rotation and reflection symmetry on the lattice, he squares and then takes a square 
root, to ensure scalar, rather than pseudoscalar, properties in the continuum limit, as in [18]. A 
torsion constraint is also necessary here. As in formulations discussed earlier, these features make 
calculations very complicated. 

Finally the authors of Ref. [14] have presented a unified treatment of Poincare, de Sitter and 
conformal gravity on the lattice. This shares many features with the formulations already described, 
so we will not discuss it further here. The main difference is that the lattice vierbein field is defined 
on the lattice links rather than at the vertices. The formulation is reflection positive, but the mode 
doubling problem seems to persist, as seen form the expansion about a flat background. 

Causal dynamical triangulations [34] are based on the action of Regge Calculus, but the ap- 
proach differs in that all simplices have identical spacelike edges and identical timelike edges, and 




(75) 



n 



^/WA]£^[il> 



Tr[ 75 P^(n) P a p(n)] exp {i\ a ^ Tr[ 7 a (l + j 5 )P^(n)}} , 

(76) 
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the discrete path integral involves summing over triangulations. In this case it is not clear how 
to use the methods discussed here and in [3], which are based on the invariant Haar measure for 
continuous rotation matrices, since this formulation does not contain explicitly continuous degrees 
of freedom which could be used for such purpose. 

The proposed formulation of Weingarten [35], based on squares, cubes and hypercubes, rather 
than simplices, involves six-index complex variables corresponding to cubes, so although it is possi- 
ble to define a large planar loop, it is not clear how to evaluate a Wilson loop, except in the special 
case when the parameter p (the coefficient of the term in the action which gives the contribution 
from the boundaries of the 4-cells) is set equal to zero, which seems to correspond to the unphysical 
case of infinite cosmological constant. 

A more radical approach to discrete quantum gravity, in which the ingredients are a set of 
points and the causal ordering between them, is known as causal sets. Recent progress includes a 
calculation of particle propagators from discrete path integrals [36]. In this formulation, it is not 
clear how to define a (closed) Wilson loop connecting points which are not causally related, and 
defining a near planar loop is also a problem here. 

6 Effects of Scalar Matter Fields 

In the next four sections, we consider whether the presence of matter affects the area law behavior 
of gravitational Wilson loops in the strong coupling limit. For each type of matter, we first describe 
briefly its transcription to the lattice [4]. 

A scalar field can be introduced as the simplest type of dynamical matter that can be coupled 
invariantly to gravity. In the continuum the scalar action for a single component field (j){x) is 
usually written as 

I[g,4>] = \ J dxy/gW d^d v <t>+(m 2 +tR)4?] + ... (77) 

where the dots denote scalar self-interaction terms. Thus for example a scalar field potential 
U {<p) could be added containing quartic field terms, whose effects could then be of interest in the 
context of cosmological models where spontaneously broken symmetries play an important role. 
The dimensionless coupling £ is arbitrary; two special cases are the minimal (£ = 0) and the 
conformal (£ = |) coupling case. In the following we shall mostly consider the case £ = 0. It 
is straightforward to extend the treatment to the case of an iV s -component scalar field (j) a with 
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One way to proceed is to introduce a lattice scalar fa defined at the vertices of the simplices. 
The corresponding lattice action can then be obtained through a procedure by which the original 
continuum metric is replaced by the induced lattice metric. Within each n-simplex one defines a 
metric 

9ij{s) = ei • ej , (78) 

with 1 < i,j < n, and which in the Euclidean case is positive definite. In components one has 
9ij = Vab e i e> j- I n terms of the edge lengths Uj = |ej — ej\, the metric is given by 

9iM = \ {ll + ll-ll) ■ (79) 

The volume of a general n-simplex is then given by 

V n (s) = ^det gij (s) . (80) 

To construct the lattice action for the scalar field, one then performs the replacement 

9i*,(x) — ► 9ij{s) 
dn<l>d u <l> — ► A^Aj<p (81) 

with the scalar field derivatives replaced by finite differences 

3^ — ► (A^) 8 = fa+^ - fa . (82) 

where the index \i labels the possible directions in which one can move away from a vertex within a 
given simplex. After some re-arrangements one finds a lattice expression for the action of a massless 
scalar field [37,38] 

/(0) = JE^(^) a . ( 83 ) 

<ij> li i 



Here vfj is the dual (Voronoi) volume [39] associated with the edge ij, and the sum is over all 
links on the lattice. Other choices for the lattice subdivision will lead to a similar formula for the 
lattice action, with the Voronoi dual volumes replaced by their appropriate counterparts for the 
new lattice. Mass and curvature terms such as the ones appearing in Eq. (77) can be added to the 
action, so that a more general lattice action is of the form 

M (fa- 4>j\ 2 , iwM/^ , tD . U 2 

<ij> 



i = \T. v!f ( M^) 2 + * £ rt d) (- 2 + ( 84 ) 
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where the term containing the discrete analog of the scalar curvature involves 

^ ( ^^E^f 2) ~V^^ (85) 

hDi 

In the expression for the scalar action, is the (dual) volume associated with the site i, and 5^ 
the deficit angle on the hinge h. The lattice scalar action contains a mass parameter m, which has 
to be tuned to zero in lattice units to achieve the lattice continuum limit for scalar correlations. 

When considering whether the gravitational Wilson loop area law holds for large loops in the 
strong coupling limit, the matter considered must be almost massless, otherwise its effects will not 
propagate over large distances and so cannot change the large Wilson loop result found in the pure 
gravity case. In fact, since the lattice Lagrangian for the scalar matter involves only factors related 
to the lattice metric (functions of the edge lengths) and not the connection (provided the parameter 
£ = 0), then the integration over the connections, which is what gives the area law, is unaffected. 



7 Effects of Lattice Fermions 

On a simplicial manifold spinor fields ifj s and tjj s are most naturally placed at the center of each d- 
simplex s. In the following we will restrict our discussion for simplicity to the four-dimensional case, 
and largely follow the original discussion given in [40,41]. As in the continuum, the construction 
of a suitable lattice action requires the introduction of the Lorentz group and its associated tetrad 
fields e^(s) within each simplex labeled by s. Within each simplex one can choose a representation 
of the Dirac gamma matrices, denoted here by 7 M (s), such that in the local coordinate basis 

ms),Y(s)} = 2sT{8). (86) 

These in turn are related to the ordinary Dirac gamma matrices 7 a , which obey 

{ 7 a , 7 6 } = 2rf> , (87) 

with rf h the flat metric, by 

*f(s) = eZ(s) 7 a , (88) 
so that within each simplex the tetrads e^(s) satisfy the usual relation 

e£(a) el(s) rf h = g^(s) . (89) 
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In general the tetrads are not fixed uniquely within a simplex, being invariant under local Lorentz 
transformations. In the following it will be preferable to discuss the Euclidean case, for which 

Vab = <W 

In the continuum the action for a massless spinor field is given by 

I = J dx^^(x)^D^(x) (90) 

where = + \ ^^ab °' ab is the spinorial covariant derivative containing the spin connection 
Lo^ab- In the absence of torsion, one can use a matrix U(V , s) to describe the parallel transport of 
any vector ^ from simplex s to a neighboring simplex s', 

P{ 8 ') = U^s',s)<j>»(s) (91) 

U therefore describes a lattice version of the connection. Indeed in the continuum such a rotation 
would be described by the matrix 

= (e Tdx Y (92) 

with rj^, the affine connection. The coordinate increment dx is interpreted as joining the center 
of s to the center of s', thereby intersecting the face f(s,s'). On the other hand, in terms of the 
Lorentz frames S(s) and S(s') defined within the two adjacent simplices, the rotation matrix is 
given instead by 

U a b (s',s) = e\{s')e\{s)U» v {s',s) (93) 

(this last matrix reduces to the identity if the two orthonormal bases E(s) and S(s') are chosen to 
be the same, in which case the connection is simply given by U (s 1 , s)^ = e^ 1 e v a ). Note that it is 
possible to choose coordinates so that U(s, s') is the unit matrix for one pair of simplices, but it 
will not then be unity for all other pairs in the presence of curvature. 

One important new ingredient is the need to introduce lattice spin rotations. Given in d di- 
mensions the above rotation matrix U(s',s), the spin connection S(s,s') between two neighboring 
simplices s and s' is defined as follows. Consider S to be an element of the 2^-dimensional repre- 
sentation of the covering group of SO(d), Spin(d), with d = 2u or d = 2u + 1, and for which S is 
a matrix of dimension 2" x 2". Then U can be written in general as 

U = exp[i(7 Q ^] (94) 

where 6 a p is an antisymmetric matrix. The cr's are \d{d — 1) d x d matrices, generators of the 
Lorentz group (SO(d) in the Euclidean case, and SO(d—l, 1) in the Lorentzian case), whose explicit 
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form is 



so that, for example, 



WapVs = ~ & p VaS 



013 



I 1 0\ 



-1 

V o ooo/ 



(95) 



(96) 



For fermions the corresponding spin rotation matrix is then obtained from 



exp i 7 a %3 



(97) 



with generators 7 a/3 = ^[7 a ,7 /3 ]- Taking appropriate traces, one can obtain a direct relationship 
between the original rotation matrix U(s, s') and the corresponding spin rotation matrix S(s, s') 



C/ a/3 = Tr(st 7a S 7/3 )/Trl 



(98) 



which determines the spin rotation matrix up to a sign. Now, if one assigns two spinors in two 
different contiguous simplices s\ and S2, one cannot in general assume that the tetrads e%(s\) and 
e^(s2) in the two simplices coincide. They will in fact be related by a matrix U(s2, s{) such that 



(99) 



and whose spinorial representation S is given in Eq. (98). Such a matrix S(s2, «i) is now needed to 
additionally parallel transport the spinor ip from a simplex s\ to the neighboring simplex S2- The 
invariant lattice action for a massless spinor takes therefore the form 



1 = \ E V (f( s ' s ')) ^ S ( U ( S > r(s') n^s, s') ^ s 

faces f(ss') 



(100) 



where the sum extends over all interfaces f(s, s') connecting one simplex s to a neighboring simplex 
s' , n^(s, s') is the unit normal to f(s, s') and V(f(s, s')) its volume. The above spinorial action can 
be considered closely analogous to the lattice Fermion action proposed originally by Wilson [1] for 
non-Abelian gauge theories. It is possible that it still suffers from the fermion doubling problem, 
although the situation is less clear for a dynamical lattice [42]. 

It is clear that the situation with gravitational Wilson loops is a bit more complicated than in 
the scalar field case, since the action now contains the spin connection matrix, which is a function 
of the matrices U which play the role of the connection. What is more, the generators of the spin 
rotation matrices are in a different representation from the generators of the rotation matrices, and 
it seems impossible to obtain, to lowest order, a spin zero object out of the combination of two 
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objects of spin one-half (S) and spin one (U), unless one applies the fermion contribution twice to 
each link, in which case a non-zero contribution can arise. We note here that if the Wilson loop 
were to contain a perimeter contribution, it would be of the form 

W{C) ~ const. {k m ) L(C) ~ exp[-TO p L(C)] (101) 

where L(C) is the length of the perimeter of the near-planar loop C, m p the particle's mass, equal 
here to m p = | In k m \ for small k m , with k m the weight of the single link contribution from the matter 
particle (sometimes referred to as the hopping parameter). Area and perimeter contributions to 
the near-planar Wilson loop would then become comparable only for exceedingly small particle 
masses, mp ~ L(C)/£ 2 , i.e. for Compton wavelengths comparable to a macroscopic loop size 
(taking A(C) « L(C) 2 /4tt). 

To demonstrate the perimeter behavior (see Fig. 10), one would need to show that the matrix 
S on the face between simplices s and s' would have a term proportional to the corresponding 
U(s,s'), with coefficient composed of 7-matrices, thereby possibly giving a non-zero contribution 
to the ^-integration. (This does not seem to be true in the infinitesimal case to lowest order, where, 
for example, S(9 34: ) = I 4 + 574 [^13 (#34) - ^24(034)]-) 




Figure 10. Illustration on how a perimeter contribution to the gravitational Wilson loop arises from matter 
field contributions. Note that now the arrows representing rotation matrices reside in principle in different 
representations . 



8 Effects of Gauge Fields 

In the continuum a locally gauge invariant action coupling an SU(N) gauge field to gravity is 

We = ~ j txyfg^fF^Fl, (102) 
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with F* u = V^Al - V v A a ^ + gf abc A b ^A c u and a, b, c = 1, . . . , N 2 - 1. On the lattice one can follow a 
procedure analogous to Wilson's construction on a hypercubic lattice, with the main difference that 
the lattice is now possibly simplicial. Given a link ij on the lattice one assigns group elements Uij, 
with each U an N x TV unitary matrix with determinant equal to one, and such that Uji = U^ 1 . 
Then with each triangle (plaquette) A, labeled by the three vertices ijk, one associates a product 
of three U matrices, 

Ua = Uijk =UijUjkUki (103) 
The discrete action is then given by [37] 

J gauge = E ^Re [Tr(l - U A )} (104) 

with 1 the unit matrix, Va the 4- volume associated with the plaquette A, A a the area of the 
triangle (plaquette) A, and c a numerical constant of order one. If one denotes by ta = cVa/Aa 
the d — 2- volume of the dual to the plaquette A, then the quantity 

A~a = °Al (105) 
is simply the ratio of this dual volume to the plaquettes area. The edge lengths lij and therefore 
the metric enter the lattice gauge field action through these volumes and areas. One important 
property of the gauge lattice action of Eq. (104) is its local invariance under gauge rotations gi 
defined at the lattice vertices., One can further show that the discrete action of Eq. (104) goes over 
in the lattice continuum limit to the correct Yang-Mills action for manifolds that are smooth and 
close to flat. 

Regarding the effects of gauge fields on the gravitational Wilson loop one can make the following 
observation. Since the action here contains no factors related to the lattice connection, the Wilson 
loop area law for large gravitational loops will be unaffected by such gauge fields, and in particular 
this will be true for the photon (which in principle could have lead to important long distance 
effects, since it is massless). 



9 Effects from a Lattice Gravitino 

Supergravity in four dimensions naturally contains a spin-3/2 gravitino, the super symmetric partner 
of the graviton. In the case of M = 1 supergravity these are the only two degrees of freedom present. 
The action contains, beside the Einstein-Hilbert action for the gravitational degrees of freedom, 



30 



the Rarita-Schwinger action for the gravitino, as well as a number of additional terms (and fields) 
required to make the action manifestly super symmetric off-shell. 

A spin-3/2 Majorana fermion in four dimensions corresponds to self-conjugate Dirac spinors 
if) n, where the Lorentz index fx = 1 . . . 4. In flat space the action for such a field is given by the 
Rarita-Schwinger term 

Crs = ~ \ e a ^ S Y£ C 75 7/3 3 7 ^ (106) 

where C is the charge conjugation matrix. Locally the action is invariant under the gauge trans- 
formation 

^V(x) ^(x) + V(x) ( 107 ) 

where e(x) is an arbitrary local Majorana spinor. 

The construction of a suitable lattice action for the spin-3/2 particle proceeds in a way that is 
rather similar to what one does in the spin-1/2 case. On a simplicial manifold the Rarita-Schwinger 
spinor fields ip p {s) and Vv(s) are mos t naturally placed at the center of each ci-simplex s. Like 
the spin-1/2 case, the construction of a suitable lattice action requires the introduction of the 
Lorentz group and its associated vierbein fields e"(s) within each simplex labeled by s, together 
with representations of the Dirac gamma matrices (see the previous discussion of Dirac fields). 

Now in the presence of gravity the continuum action for a massless spin-3/2 field is given by 

h/2 = ~\ J dx^ge^ Xa ^^x) lblv D x ^ a {x) (108) 

with the Rarita-Schwinger field subject to the Majorana constraint W = Ctp^x) 7 '. Here the 
covariant derivative is defined as 

D^p = 0^ p - Tl p ^ + \ uj vah a ab W (109) 

and involves both the standard affine connection as well as the vierbein connection 

w va b= \ [e/idv e bfi ~ e bu ) + e a p e b a {d a e cp ) e c u ] 

- (a++b) (110) 

with Dirac spin matrices a a b = ^[j a , 7b], and e^ 90 the usual Levi-Civita tensor, such that e^ pcr = 

_ gel xvpa 

It is easiest to just consider two neighboring simplices s\ and S2, covered by a common coordinate 
system x^. When the two vierbeins in si and «2 are made to coincide, one can then use a common 
set of gamma matrices 7 M within both simplices. Then (in the absence of torsion) the covariant 
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derivative in Eq. (108) reduces to just an ordinary derivative. The fermion field Vv( x ) within 
the two simplices can then be suitably interpolated, and one obtains a lattice action expression 
very similar to the spinor case. One can then relax the condition that the vierbeins e^(si) and 
e^(s2) in the two simplices coincide. If they do not, then they will be related by a matrix U(s2, si) 
such that 

<{S2) = U» u {s 2 ,s l )e» a {s l ) (111) 

and whose spinorial representation S was given previously in Eq. (98). But the new ingredient in 
the spin-3/2 case is that, besides requiring a spin rotation matrix S(s2, si), now one also needs the 
matrix U^(s,s') describing the corresponding parallel transport of the Lorentz vector Vv( s ) from 
a simplex s\ to the neighboring simplex S2- The invariant lattice action for a massless spin-3/2 
particle takes therefore the form 

1 = -\ E ^(/( S ,/))e^^( S )S(U( S , S '))7,(s')nA( S ,^^ CT ^« / )V'p(s / ) (H2) 

faces f(ss') 

with 

M S )S(XJ( S ,s')) 7v (s')Ms') = ^a(8)S a p (XJ(8,s')) lv ^(s')^ p (8') (113) 
and the sum Xfaces f(ss') extends over all interfaces f(s, s') connecting one simplex s to a neighboring 
simplex s'. When compared to the spin-1/2 case, the most important modification is the second 
rotation matrix U u fl (s,s'), which describes the parallel transport of the fermionic vector ip^ from 
the site s to the site s', which is required in order to obtain locally a Lorentz scalar contribution 
to the action. 

In this case again one expects the Wilson loop to follow a perimeter law, as in the spin one-half 
case of Eq. (101), because the matter action explicitly contains factors of U which will contribute 
when the Us and Ss around the loop are integrated over, which of course requires that one also 
take into account the spin connection matrices. These add complexity but are not expected, due to 
the nature of the interaction, to change the answer. The same general considerations then apply as 
in the spin-1/2 case: the perimeter contribution to the gravitational Wilson loop can significantly 
modify the area law result only if the corresponding particle mass is exceedingly small. 

10 Physical Consequences 

Having ascertained with some degree of confidence that in a number of different, and quite un- 
related, lattice discretization of gravity the gravitational loop follows an area law for sufficiently 
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strong coupling G, 

<W(C)> ~ exp(-A c /f) (114) 

with £ determined by scaling to be the nonperturbative gravitational correlation length, let us now 
turn to possible physical interpretations of the results. In particular, we will explore their connection 
to a semi-classical picture, describing the properties of curvature on very large, macroscopic distance 
scales Note that, in the above expression, £ is intended to be the renormalization group invariant 
quantity obtained by integrating the /3-function for the Newtonian coupling G, 

CHG) = const. A exp (- f . (115) 

with A the ultraviolet cutoff (and analogous to Eq. (8) for gauge theories). Then in the vicinity of 
the ultraviolet fixed point at G c 

(3(G)) = A* G(ji) ^ P'(G C ) (G — G c ) + . . . , (116) 

which gives 

r'CGJcx A\(G-G C )/G C \" , (117) 

with a correlation length exponent v = —1//3'(G C ). In particular the correlation length £(G) is 
related to the bare Newtonian coupling G, and diverges, in units of the cutoff A, as one approaches 
the fixed point at G c . Thus for a bare G very close to G c the two scales, A and S," 1 can be vastly 
different. Furthermore the result of Eq. (114) was derived from the lattice theory of gravity in the 
strong coupling limit G — > oo. But one would expect, based on general scaling arguments, that such 
a behavior would persist through the whole strong coupling phase G > G c , all the way up to the 
non-trivial ultraviolet fixed point at G c . This is indeed what happens in non-abelian gauge theories 
and spin systems: the only scale determining the non-trivial scaling properties in the vicinity of 
the fixed point is £. The corresponding behavior is then known as universal renormalization group 
scaling. 

As discussed at the beginning of this work and in [3,4], the rotation matrix appearing in the 
gravitational Wilson loop can be related classically to a well-defined classical physical process, one 
in which a vector is parallel transported around a large loop, and at the end is compared to its 
original orientation. Then the vector's rotation is directly related to some sort of average curvature 
enclosed by the loop; the total rotation matrix U(C) is given by a path-ordered (V) exponential 
of the integral of the affine connection T^ u , as in Eq. (1). In a semi-classical description of the 
parallel transport process of a vector around a very large loop, one can re-expresses the connection 



33 



in terms of a suitable coarse-grained, semi-classical slowly varying Riemann tensor, as in Eq. (2). 
Since the rotation is small for weak curvatures, one has 

U%{C) ~ [l + \ R.^AZ + ...] a p . (118) 

At this stage one can compare the above semi-classical expression to the quantum result of Eqs. (44), 
(66) and (114), and in particular one would like to relate the coefficients of the area terms. Since 
one expression [Eq. (118)] is a matrix and the other [Eq. (114)] is a scalar, one needs to take the 
trace after first contracting the rotation matrix with (Be + eh), as in our second definition of the 
Wilson loop of Eq. (16), giving 

W(C) ~ Tr ((B c + eh) exp j \ R-.^Ag . (119) 

Next, it is advantageous to consider the lattice analog of a background classical manifold with 
constant or near-constant large scale curvature, 

RnuXa = I A (g^ u g\ a - g^x 9ua) (120) 

so that here one can set for the curvature tensor 

Ra (3nv = RB a /3 B fMU , (121) 

where R is some average curvature over the loop, and the area bivectors B here will be taken 
to coincide with Be- The trace of the product of (Be + eh) with this expression then gives 
Tr(R B c Ac) = — 2 R Ac- This is to be compared with the linear term from the other 
exponential expression, — Ac/£ 2 - Thus the average curvature is computed to be of the order 

R ~ +l/£ 2 (122) 

at least in the small k = 1/8itG limit. An equivalent way of phrasing the last result makes use of 
the classical field equations in the absence of matter R = 4A. Then l/£ 2 should be identified, up 
to a constant of proportionality of order one, with the observed scaled cosmological constant A (, s , 

Kbs - + 4 ■ ( 123 ) 

The latter can be regarded either as a measure of the vacuum energy, or of the intrinsic curvature of 
the vacuum. One sees therefore that a direct calculation of the gravitational Wilson loop provides a 
direct insight into whether the manifold is De Sitter or anti-De Sitter at large distances. Moreover, 
in the case of lattice gravity at strong coupling, as has been shown in this work, it seems virtually 
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impossible to obtain a negative sign in Eqs. (122) or (123), which would then suggest that quantum 
gravity predicts a positive cosmological constant at large distances (Again, the analogy with non- 
abelian gauge theories comes to mind, where one has for the non-perturbative gluon condensate 
< > ~ I/? 4 ) where £ is the non-perturbative QCD correlation length, £qcd ~ A^; the 
analog of the vacuum condensate in non-abelian field gauge theories is then naturally seen here to 
be the vacuum expectation value of the curvature) . 

At first it would seem that in principle the scale £ could take any value, including very small ones, 
based on the naive estimate £ ~ lp, where lp is the Planck length whose magnitude is comparable 
to the (inverse of the) ultraviolet cutoff A. This of course would then preclude any observable 
quantum effects in the foreseeable future. But the relationship between £ and large scale curvature, 
or more precisely between £ and A ;, s , arising out of the specific properties of the gravitational 
Wilson loop, would seem to imply a strikingly large, cosmological value for £ ~ 10 28 cm, given the 
present observational bounds on X Q b s . Closely related possibilities exist, such as an identification 
of £ with the Hubble constant as measured today, £ ~ 1/Hq; since this quantity is presumably 
time-dependent, a possible scenario is one in which = = lim^oo H(t), with = | A f> s . 

One should note here that the link between the long distance cosmological constant \ t, s and the 
renormalization group invariant mass scale m = l/£ is not entirely unexpected, and can be viewed 
from a different perspective. If one compares for example, at the most basic level, the structure of 
the three classical field equations for gravity, QED (made massive via the Higgs mechanism) and 
a self-interacting scalar field, respectively, 

Rfiu — \ g^u R + A g^y = 8nG 
dT^ + m 2 A u = 4vre j u 

d^d^ + m 2 ^ = !</> 3 (124) 

one finds that it is indeed consistent and natural to regard the (scaled) cosmological constant as 
a renormalization group invariant quantity, just like the particle's masses in the other two field 
theories. It should also be stressed that in all cases, the scale £ or m is not predicted from first 
principles, instead its value must be taken from observation and experiment. In the case of gravity, 
all one can say is that a specific value for £ is related to how close the bare coupling G(A) is to the 
ultraviolet fixed point value G c , as shown for example in Eq. (117). The factor of G — G c there 
is what accounts for the possibility of having in quantum gravity two widely distinct scales, the 
cutoff A and the non-perturbative scale £. 

Finally it is important to note though that another important and independent physical pre- 
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diction arises from the identification of £ with 1/V \ bs '■ one expects £ to determine the scale 
dependence of the couplings appearing on the r.h.s. of the field equations in Eq. (124), and in 
particular the scale dependence of Newton's constant G(fi), as obtained from solving for example 
the renormalization group equations for G, Eqs. (115) and (116). 

As discussed in [44], a running of the gravitational constant of the type discussed in [7] is 
phrased best in a fully covariant formulation, such as an effective classical, but non-local, field 
equation of the type 

iV -\o^R + \g^ v = totG{u)Tp, (125) 

with A ~ l/£ 2 , and G(U) the running Newton's constant 

G -► G(p) (126) 

with the running given by 



<?(□) = G c 



1 \ 2i/ 



I + a ( ^75 ) + 



(127) 



and ao ^ 42 > and v ~ 1/3 [45]. G c in the above expression should be identified to a first 
approximation with the laboratory scale value \[G~ C ~ 1.6 x 10 -33 cm [44,4]. The running of G can 
then be worked out in detail for specific coordinate choices. Thus in the static isotropic case one 
finds a gradual slow increase in G with distance, in accordance with the formula 



G -► G(r) = G 1 + — m 3 r 3 In -— ^ + . . . (128) 



a o s s i 1 
— m° r° In — — 

3 7r m z r z 

in the regime r»2MG, where 2MG is the horizon radius, and m = l/£. This would then suggest 
that the natural scale entering the quantum scale dependence of G(r) is the observed cosmological 
constant, £ = l/V A D & S . 



11 Conclusions 

From our study of Wilson loops, where defined and calculable, in all theories of discrete gravity 
that we have found, it seems that the area law holds for large loops in the strong coupling domain. 
This suggests that we can infer, as in [3], that a universal prediction of strongly coupled gravity 
without matter is that the scaled cosmological constant is positive. 

The effect of adding matter coupled to gravity is less clear-cut, although it is only massless or 
almost massless matter which propagates to sufficiently large distances to affect large gravitational 
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Wilson loops. In that case, scalar matter and gauge fields (in particular the photon) do not affect 
the area law. For very low mass fermions (e.g. neutrinos), it is possible that the coupling gives 
rise to a perimeter contribution, which could replace the area law for suitable ratios of coupling 
constants, but this seems unlikely. Similarly, the lattice gravitino could produce a perimeter law. 
These possibilities will be investigated in future work. Numerical simulations of simplicial lattice 
gravity could provide vital clues here [45]. 
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